Second harmonic microscopy of monolayer M0S2 
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We show that the lack of inversion symmetry in monolayer M0S2 allows strong optical second 
harmonic generation. Second harmonic of an 810-nm pulse is generated in a mechanically exfoliated 
monolayer, with a nonlinear susceptibility on the order of 10~ 7 m/V. The susceptibility reduces 
by a factor of seven in trilayers, and by about two orders of magnitude in even layers. A proof- 
of-principle second harmonic microscopy measurement is performed on samples grown by chemical 
vapor deposition, which illustrates potential applications of this effect in fast and non-invasive 
detection of crystalline orientation, thickness uniformity, layer stacking, and single-crystal domain 
size of atomically thin films of M0S2 and similar materials. 



Recently, there is a growing interest in exploring new 
types of atomically thin crystals based on layered ma- 
terials, such as transition metal dichalcogenides.— The 
most extensively studied member of this family is M0S2. 
In 2010, photoluminescence experiments 2 ^ and micro- 
scopic calculations^ indicated that, although bulk M0S2 
is an indirect semiconductors, its monolayer is a direct 
semiconductor with a bandgap of about 1.88 eV. Such a 
large bandgap and the structural similarity with widely 
studied graphenc immediately stimulated investigations 
on its potential applications in logic electronics.^ In 
2011, top-gated transistors based on M0S2 monolay- 
ers were fabricated^ Later, integrated circuits based on 
monolayer— and bilayer— for logic operations were demon- 
strated. Ambipolar transport in a multilayer transistor 
gated by ionic liquids was also demonstrated, showing the 
feasibility to develop p-n-junction-based devices* 1 ^ Fur- 
thermore, the recently demonstrated superior strength 
and flexibility 1 ^ make M0S2 atomic layers an attrac- 
tive candidate for applications in flexible electronicsi 12 ' 13 
In addition, since monolayer M0S2 has a bandgap in 
the visible range, has workfunctions that are compatible 
with commonly used electrode materials, and has sta- 
ble charge exciton state even at room temperature,— it 
is also an attractive candidate for various optoelectronic 
and photonic applications, such as phototransistor a 15 ' 16 
and hctcrojunction solar cells.— 

In contrast to these breakthroughs in understanding 
the mechanical, electronic, and linear optical properties 
of M0S2 atomic layers, little is known about their non- 
linear optical properties. Nonlinear optical responses are 
important aspects of light-matter interaction, and can 
play important roles in various photonic and optoelec- 
tronic applications, especially in those involving high in- 
tensity laser beams. Bulk M0S2 crystal with 2H stack- 
ing order belongs to space group Dq^, which is inver- 
sion symmetric. Hence, its second-order nonlinear re- 
sponse should vanish.— Indeed, one early experiment 
showed that second-order nonlinear susceptibility of 2H 
bulk M0S2 is at most 10~ 14 m/V.— However, the in- 
version symmetry is broken in a monolayer, which has 
D%h symmetry. One consequence of such a symme- 
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FIG. 1. (a) Schematics of the experimental setup, (b) Mi- 
croscope image of a mechanically exfoliated M0S2 flake, (c) 
Lattice structure of monolayer M0S2 • 



try reduction is to allow valley- selective optical inter- 
band transitions, which has been observed by several 
groups recentl y 14 ' 20 " — and can be used for valley tronics, 
in which the valley index of electrons is used to carry in- 
formation. Here we show that the lack of inversion sym- 
metry allows unusually strong optical second harmonic 
generation (SHG) in monolayer M0S2 flakes prepared 
by mechanical exfoliation and chemical vapor deposi- 
tion (CVD). This effect is very sensitive to layer thick- 
ness, crystalline orientation, and layer stacking. Based 
on these properties, we demonstrate a second harmonic 
microscopy for characterization of M0S2 thin films. 

Figure nja) shows the experimental setup. The fun- 
damental pulse with an angular frequency uj and a cen- 
tral wavelength of 810 nm is obtained from a Tksapphirc 
laser. It is tightly focused to a spot of 2 /im (full width 
at half maximum) by a microscope objective lens. The 
second harmonic (SH) generated is collected by the same 
lens, and detected by a spectrometer equipped with a 
thermoelectric cooled silicon charge-coupled device cam- 
era. A set of color filters is used to block the fundamen- 
tal and other unwanted light. With 810-nm wavelength, 
band-to-band absorption of fundamental is avoided, and 
both fundamental and SH can be detected efficiently with 
silicon detectors, which facilitates alignment and location 
of the laser spot. Figure [ljb) is a microscope photo of a 
M0S2 flake that is mechanically exfoliated onto a Si/SiC>2 
(90 nm) substrate. The region marked with red label 
1L is identified as a monolayer according to its optical 
contrast^ — Raman spectrum ; 28 ' 29 and photolumincs- 
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ence spectrumj2i£ Other regions with few atomic layers, 
as indicated by the red labels, are assigned according to 
their relative optical contrasts. 

The structure of monolayer M0S2 is schematically 
shown in Fig. [TJc), where each yellow circle represents 
two S atoms vertically separated by 0.65 nm, and blue cir- 
cles indicate plane of Mo atoms located between the two S 
atomic planes. With the D^h symmetry, the second-order 
nonlinear susceptibility tensor has nonzero elements of 
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armchair direction, which is 30° from the zig-zag direc- 
tion, along which the mirror symmetry is broken. In 
the experiment, the fundamental beam is normal inci- 
dent (along —z') and is linearly polarized along horizon- 
tal direction [defined as x in the laboratory coordinates, 
as shown in Fig. [T{b)]. It is straightforward to show 
that the parallel (x) and perpendicular (y) components 
of SH field are proportional to sin36> and cos3#, respec- 
tively, where 9 is the angle between x and x' . Hence, the 
power of the two components varies as P x oc sin 2 3f? and 
P y oc cos 2 3#, while the total power is independent of 9. 

In our experiment, we first measure the total SH power 
with a fundamental power of 4 mW. The upper inset of 
Fig. [2ja) shows the spectra of the SH (blue) and the fun- 
damental (red), confirming that the former is indeed at 
half wavelength of the latter. The gray curve is a spec- 
trum (multiplied by a factor of 100) taken under the same 
conditions but with the laser spot located on bare sub- 
strate. Hence, the contrast of monolayer with respect to 
substrate is at least 10 4 , which is much higher than lin- 
ear optical microscopy (about 0.3). The main panel of 
Fig. [Ha) shows how the SH power varies with the fun- 
damental power. The peak irradiance of fundamental 
and SH pulses, deduced from the powers, are also plot- 
ted for convenience, as top and right axes. The solid 
line is the expected quadratic dependence for the SHG 
process. Next, by placing a linear polarizer in front of 
the spectrometer, we measure P x and P y as a function 
of 9, the angle between x and x' , by rotating the sample 
about z axis. Figure [2jb) shows the results, along with 
the expected 9 dependence (solid lines) from the D 3H 
symmetry. 

In order to estimate the magnitude of x*- 2 ' from the 
measurement, we model the monolayer as a bulk medium. 
Since the flake thickness (d = 0.65 nm) is much smaller 
than the coherence length, the SHG is not influenced 
by phase-matching conditions. By solving the coupled- 
wave equations,— the SH field amplitude of the parallel 
component 
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where c is the speed of light in a vacuum and n% u is 
the index of refraction at SH, and £ u is the fundamen- 
tal field amplitude. The £ x is related to the irradiance 
by I x = n2 UJ eoc£ x £*/2, which can be calculated from the 
measured quantity, average power, by considering that 
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FIG. 2. Second harmonic generation from mechanically exfo- 
liated M0S2 sample: The upper inset of (a) shows the spectra 
of second harmonic from the monolayer M0S2 and from bare 
substrate (gray, multiplied by a factor of 100), as well as the 
fundamental beams (red). The lower inset shows the second 
harmonic power measured from regions with different atomic 
layers. The main panel of (a) shows the power dependence 
of second harmonic generation, with the solid line indicat- 
ing the expected quadratic dependence, (b) Power of parallel 
(blue squares) and perpendicular (black circles) components 
of second harmonic as a function of 6, the angle between the 
laboratory and the crystalline coordinates. The blue (black) 
solid line indicates the expected sin 2 3f5 (cos 2 3f?) dependence. 



I x is Gaussian in both time and space, with widths (full 
width at half maxima) of r and W, respectively By using 
W = 2 fim, t = 200 fs, / = 81 MHz, n 2uj « 6.0,-21 an d re- 
flection coefficient of 0.09 from this multilayer structure, 
we find that the magnitude of x 1S about 10 -7 m/V. 
We note that due to the nonlinear nature of this process, 
such a deduction replies on accurate knowledge on many 
experimental parameters, such as the shape and dura- 
tion of the fundamental pulse, the shape and size of the 
focused fundamental spot at sample, and the relation be- 
tween the measured spectral counts and the actually SH 
power. Hence, this value should be viewed as an order- 
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of-magnitudc estimate. However, the relative comparison 
of throughout this paper are not influenced by such 
uncertainties, and are thus accurate. 

Since monolayer M0S2 possesses such a large x^ 2 \ 
which vanishes in bulk, it is interesting to study how 
t^( 2 ) varies with the number of atomic layers. We mea- 
sure the total power of SH from different regions of the 
flake shown in Fig.[TJb), with a fixed fundamental power 
of 4 mW. The results are summaries in the lower inset 
of Fig. EJa). Since the total power is independent of 6, 
the measurement is not influenced by potentially different 
crystal orientations of these regions. We find that x^ of 
trilaycr is about a factor of seven smaller than monolayer, 
while those of bilayer and quadralayer are about two or- 
ders of magnitude smaller than the monolayer. Since 
flakes with even number of atomic layers possess inver- 
sion symmetry, their second-order response should van- 
ish. The smaller but nonzero x^ can be attributed to 
surface and interface effects. According to this measure- 
ment, the contrast of monolayer with respect to bilayer 
and quaralayer is about 10 4 . Similar layer- number de- 
pendence has also been observed recently in WS2 and 
WSe2 £2- We also measure a thick flake that can be con- 
sidered as a bulk. The SH power is about five orders of 
magnitude smaller than the monolayer, indicating a very 
small x^ 2 \ as a result of the inversion symmetry. 

The observed SHG can be used to fast and in-situ char- 
acterize atomically thin films of M0S2 and similar mate- 
rials. Although high quality monolayer M0S2 can be pro- 
duced by simple mechanical exfoliation^ 1 - and identified 
by optical contrast with certain substrate a 25 ' 26 and Ra- 
man spectroscopy ) 28 ' 29 applications of this material rely 
on development of scalable techniques. Following initial 
works of mechanical exfoliation, other top-down methods 
with better potential for large-scale production have been 
developed, such as lithium ion exfoliation^— and ul- 
trasonic exfoliation in liquidsi^Zr— Promising progresses 
have also been made in developing bottom-up methods, 
including hydrothermal synthesi a 44 ' 45 and CVD on insu- 
lating substrates 4 ^— and graphene^ However, one sig- 
nificant obstacle is the lack of techniques for fast and in- 
situ sample characterization. For example, thin films of 
M0S2 fabricated by these techniques are polycrystallinc. 
They are composed of single-crystalline domains with 
random crystal orientations and are separated by grain 
boundaries, which severely limit performance of the films, 
especially their conductivity and mechanical strength. 
However, it is difficult to locate the grain boundaries and 
monitor size of these domains in-situ. 

The lower inset of Fig. [2ta) shows a microscope 
photo of some triangular monolayer M0S2 flakes on a 
Si/SiO 2 (280 nm) substrate fabricated by CVD. The sam- 
ples were prepared using M0O3 and sublimated sulfur as 
precursors. M0O3 is positioned close to the designated 
growth substrate at the center of the furnace, while sub- 
limated sulfur is positioned upstream at a zone where 
evaporation starts at 750°C. The reaction of the precur- 
sors at 850°C in a furnace flushed with nitrogen results in 
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FIG. 3. Second harmonic generation from a triangular mono- 
layer M0S2 flake grown by CVD, as shown in the lower inset 
of (a). The main panel of (a) shows the power dependence of 
second harmonic generation. The solid line indicates the ex- 
pected quadratic dependence. Panel (b) shows angular depen- 
dence of the parallel (blue squares) and perpendicular (black 
circles) components of second harmonic, along with the ex- 
pected dependence (solid lines). The upper inset of (a) shows 
a separate measurement of the parallel component with a finer 
step size near = 0°. 



nucleation of single crystalline domains. The density of 
nucleation and samples sizes can be controlled by mon- 
itoring the pressure and the closely related sulfur con- 
centration in the chamber. By maintaining a positive 
pressure in the range of 5 - 20 KPa, M0S2 domains with 
sizes in the range of 10 - 40 [im. are synthesized, with a 
ramping time of 60 to 90 minutes and 10 minutes at the 
reaction temperature^ 1 - 

The main panel of Fig. (3{a) shows the quadratic 
power dependence of SHG, similar to Fig.[^a), measured 
from the well-separated flake on which the crystalline 
and laboratory coordinates (x' and x, respectively) are 
plotted [the lower inset of Fig. [2ta)]. By rotating the 
sample, we measure P x and P y as a function of 9, as 
shown in Fig. [3Jb) . The results are similar to the exfo- 



liated sample shown in Fig. [2jb). A separate measure- 
ment of the parallel component with finer resolution near 
9 = 0°, shown in the upper inset of Fig. EJa), confirms 
that the minimal parallel component occurs precisely at 
8 = 0°. The edges of these triangular flakes are expected 
to be along zig-zag directions since these are lowest en- 
ergy configurations . 52 ' 53 The maximum parallel compo- 
nent of SH should occurs when the fundamental is polar- 
ized along the zig-zag direction, which is consistent our 
observation. Hence, the SHG further confirms that the 
direction of the edges is zig-zag. We repeat the measure- 
ment with several other similar triangular flakes, and ob- 
tained the same result. Such an established relation also 
allows us to determine the crystal orientation of the me- 
chanically exfoliated sample shown in Fig. [ljb): that is, 
the armchair direction of the 1L region is horizontal, and 
its lattice orientation is as shown in Fig. [TJc). From the 
strength of the SH, we deduce a x' 2 ' w 5 x 10~ 9 m/V, 
which is about a factor of 20 smaller than the mechani- 
cally exfoliated flake. 

Figure [4] summarizes our proof-of-principle demonstra- 
tion of a polarization-revolved SH microscopy. We study 
a region on the substrate with quasi-continuos films, as 
shown in Fig. HJa). It contains a high density but still 
separated and randomly oriented triangular flakes, so 
that we can correlate domains observed in SHG to the 
actual regions. In this measurement, we scan a 20-mW 
fundamental spot across the region indicated by the box 
in Fig.|4|a), and detect the powers of the parallel and per- 
pendicular components of the SH, as shown in Figs. EJb) 
and (c), respectively. The errors in these scans are be- 
low 2 pW, or smaller than 1% of the maximum signal. 
Figure 0Jd) shows the total power, obtained by adding 
(b) and (c). From (b) and (c), we calculate the angle 
by using 8 = (l/3)tan _1 y / P :E /Pj,, as shown in (e). The 
uncertainty on the angle is below 1°. 

The combination of linear and nonlinear optical mi- 
croscopy can provide valuable information on polycrys- 
talline thin films grown by CVD. First, the greenish dots 
in Fig. QJa) indicate that at the central area of some 
flakes, a second (or even third) layer is grown. Second 
harmonic images show that these areas have higher SH 
power. In bilaycr M0S2 exfoliated from 2_ff-stacked crys- 
tals, the two layers are inversely oriented so that the bi- 
layer possesses inversion symmetry. Hence, it second- 
order response should vanish, as confirmed in Fig. HJa) . 
The higher SH power observed from multilayer regions 
of CVD-grown flakes indicates that these multilayers are 
not 2i7-stacked. This is similar to multilayer graphene 
grown by CVD. Clearly, the SH microscopy is capable of 
probing relative orientations among multilayers of M0S2. 
Second, panel (e) shows that 8 is uniform over the left 
flake, which is about 15°. This is consistent with the 
shape observed in (a) (white dashed line). With further 
growth time, this flake will merge with other flakes to 
form a continuous polycrystalline film. Linear optical 
microscopy would not allow identification of each single- 
crystalline domains. However the 8 map can still distin- 
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FIG. 4. (a) Optical microscopy photo of a region of substrate 
containing flakes grown by CVD. (b) and (c): Maps of P x and 
P y over the region indicated by the box in (a), (d) Map of 
the total power, P x + P y . (e) Map of 6 calculated from (b) 
and (c). 



guish these domains, locate their boundaries, and mea- 
sure their sizes. Third, the shapes of these flakes are ir- 
regular in the SH maps. Especially, the parallel and per- 
pendicular components have different edge shapes. This 
can be attributed to the roughness on the edges and dif- 
ferent termination configurations on the edges. Although 
further characterizations are needed to correlate the mi- 
croscopic structure on the edges to the SH power, this ob- 
servation illustrates the potential of using SHG to study 
these edge structures. Finally, although the three flakes 
look similar in (a), the SH power are different, and the 
9 of the middle and right flakes are irregular. This illus- 
trates that the SH microscopy can show different prop- 
erties and qualities of the flakes that the linear optical 
microscopy cannot. However, further studies are needed 
to correlate SHG to these specific sample characteristics. 

In summary, we have observed strong second harmonic 
generation in monolayer M0S2 fabricated by mechanical 
exfoliation and CVD, and performed a proof-of-principle 
second harmonic microscopy measurement. Our results 
show that such a nonlinear optical effect can be used to 
fast and non-invasively characterize atomically thin films 
of M0S2 and other similar materials. Compared to lin- 
ear optical microscopy, the contrast is enhanced by at 
least four orders of magnitude. Unlike linear optical mi- 
croscopy that replies on light interference from carefully 
designed multilayer substrates, such a nonlinear optical 
microscopy can be applied to any substrates with weak 
second-order nonlinearity, such as silicon and glass. Al- 
though Raman microscopy has also been used to identify 
monolayer M0S2, the Raman shift often depends on sub- 
strates, and the contrast is relatively low. In addition 
to these advantages in identifying monolayers, the sec- 
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ond harmonic microscopy can probe crystal orientation, 
single-crystal domain size, and layer stacking. 
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